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Abstract: We investigate the construction of tree-level MHV gluon amplitudes in multiplet bases using 
BCFW recursion. The multiplet basis decomposition can either be obtained by decomposing results derived 
in (for example) the DDM basis or by formulating the recursion directly in the multiplet basis. We focus 
on the latter approach and show how to efficiently deal with the color structure appearing in the recursion. 
For illustration, we also explicitly calculate the four-, five- and six-gluon amplitudes. 

Keywords: QCD, Multiplet bases, Recursion relations, Scattering amplitudes 


1 


On leave from Center for Field Theory and Particle Physics, Department of Physics, Fudan University. 




Contents 


1 Introduction 1 

2 Evaluation by scalar products 6 

2.1 The three-gluon amplitude 7 

2.2 The four-gluon amplitude 7 

3 Recursion in multiplet bases 9 

3.1 BCFW recursion for color-dressed amplitudes 10 

3.2 Kinematic recursion 12 

3.3 Color structure recursion 14 

3.3.1 Example: 4 —> 5 gluons 17 

3.3.2 The general case: n — 1 —> n gluons 20 

3.4 The five-gluon amplitude 24 

3.5 The six-gluon amplitude 26 

4 Conclusion and outlook 26 

A Vanishing of the (3,n-l) division 27 

B Five-gluon multiplet basis, radiation matrices and MHV amplitudes 28 


1 Introduction 

Scattering amplitudes are essential tools for understanding multi-parton processes in high energy physics. 
Traditional calculations based on Feynman diagrams have proven extremely cumbersome due to the large 
number of diagrams and the correlation of the color and kinematic factors. To simplify the calculations 
of amplitudes, large efforts have been made, including the development of the spinor helicity formalism 
[1-5], the Berends-Giele recursion for off-shell currents [6], various forms of color decompositions [7-27], 
the Parke-Taylor formula for color-ordered maximally helicity violating (MHV) amplitudes [6, 28], the 
Kleiss-Kuijf (KK) relation [29], the twistor string method [30], Cachazo-Svrcek-Witten (CSW) rules [31], 
Britto-Cachazo-Feng-Witten (BCFW) recursion [32, 33] and the Bern-Carrasco-Johansson (BCJ) relation 
[34] between color-ordered amplitudes. Reviews of related topics can be found in [35-43]. 

In most of this work the color decomposition plays an important role for understanding and simplifying 
scattering amplitudes. Several types of color decompositions for gluon amplitudes are available, most 
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notably the trace basis decomposition [11-20], the Del Duca-Dixon-Maltoni (DDM) basis decomposition 
[21, 22] and the color-flow basis decomposition [23, 24], Since recently it is also known how to construct 
orthogonal group theory based multiplet bases for any number of quarks and gluons [25]. 

While the topic of scattering amplitude recursion relations has been explored for the kinematic factors 
(the color-ordered amplitudes) of the first three bases for a while [6, 32, 33], this field is unknown territory 
in the context of multiplet bases. In the present paper we take the first step to remedy this by showing how 
to use BCFW recursion relations for multiplet bases in the case of MHV amplitudes in pure Yang-Mills 
theory. To set the stage for this task, we first give a brief overview of the standard color decompositions 
and the present status of recursion strategies. 


• Trace bases 

The color decomposition for a tree-level gluon amplitude with n gluons in trace bases is given by 1 
[11-19], 


M(gi,g 2 , ...,g n )=g 


n —2 


i 

7 ^, 


Y Tv(t 91 t 9 ^ ...t 9 ^)A(a), 

a, s.t. 01 = 1 


( 1 . 1 ) 


where we have used a general generator normalization tr(t a t b ) = T^6 ab . The cyclicity of the trace 
allows for fixing a\ = 1, thus leaving (n— 1)! color structures in the sum. The kinematic factors, A(a) 
are called color-ordered amplitudes and can be calculated from the color-ordered Feynman rules, also 
known as the color-stripped Feynman rules. These bases (spanning sets) are easily extended to loop 
level, upon which products of traces occur, [16, 17, 19, 44], as well as to processes with quarks, 
requiring open quark-lines in addition to the traces [14-17, 19, 35]. 


• Color-flow bases 

An approach similar to the trace basis approach is given by the color-flow bases [23, 24], Here the 
gluon field is rewritten in terms of the fundamental representation ( A At )*-, i,j = 1, ..,N C , for N c colors, 
and the color structure is described in terms of flow of color. For tree-level gluon amplitudes the 
color decomposition is given by 


M(g 1 ,g 2 ,...,g n ) = g n 2 



£ 

<r, s.t. <ti = 1 


Jo - 2 




( 1 . 2 ) 


where the sum runs over the (n — 1)! permutations from connecting color lines. It is not hard to 
argue that the amplitudes A(a) are the same as in the trace bases. Similar to the trace bases, these 
bases are extendable to processes at higher order. Their advantage lies in better scaling properties 
for Monte Carlo treatment of the color structure [24]. 

1 In the color decomposition formulae, we suppress the helicity of the external legs for convenience. Only when discussing 

amplitudes with a particular helicity configuration, we specify the helicity information of the external legs. 
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• Del Duca-Dixon-Maltoni bases 


Tree-level gluon amplitudes may alternatively be decomposed using the Del Duca-Dixon-Maltoni 
(DDM) bases [21, 22] 

n —2 

_ i p n -ag° n -i9* n A (e), (1.3) 

cr, s.t. (Tp — 1, cr n =n 

where (j\ = 1 and a n = n are fixed. The color-ordered amplitudes A{a\ = 1,02, ... ,a n -\,a n = n) 
form the so called Kleiss-Kuijf (KK) basis [29]. All other color-ordered amplitudes, i.e., all the 
amplitudes in eq. (1.1) with cr(n) ^ n , can be expressed in the KK basis using the Kleiss-Kuijf (KK) 
relation 


M(gi,g 2 , ...,g n ) = g 


n —2 


i 

7^. 


A{l,{a},n,{P}) = {-l) nf> ^2 A(l, a, n), (1.4) 

crdOP(a IJ 0 T ) 

where ng denotes the number of indices in the set /3, and the sum runs over all possible permutations 
which keep the relative order of indices in the index set a and reverse the relative index order in the 
index set /3, while allowing for all possible relative orderings between each on with respect to each /3j. 

Compared to the trace bases and color-flow bases there is clearly an advantage in needing only 
(n — 2)! rather than (n — 1)! spanning vectors. The color decomposition in DDM bases can also be 
extended to one-loop level and to amplitudes containing a quark-antiquark pair [22], A proof of the 
KK relation can be found in [22], the BCFW approach to the KK relation is presented in [45], and 
the Berends-Giele recursion approach for the off-shell KK relation was given in [46]. 

In the three color decompositions above, the kinematic factors can be expressed using the color-ordered 
amplitudes that were defined in the trace bases decomposition, eq. (1.1). Employing BCFW recursion, 
one can express these amplitudes in terms of products of lower point on-shell color-ordered amplitudes. 
Specifically, for the n-gluon color-ordered tree amplitude A(l, 2, 3,..., n), if we shift the momenta of gluon 
1 and n using some complex four-vector q and some complex variable z, 

pi(z)=pi-zq, p n (z)=p n + zq, (1.5) 

such that the gluons 1 and h remain on-shell, the color-ordered amplitude A(l, 2, 3,..., n) is given by the 
BCFW recursion (see figure 1) 

A (1, 2,..., n) 

n —2 

= EE A (^)> 2 ’ 

i =2 h=± 

i 

where, -Pi, 2 ,...,i = Pk- F° r a given i, Zi can be solved from 
fc=i 


P l(zi),2,...7 Zi 


P 2 

1,2 


-A 


5 -h 


PA n _ .,* + 1 ,. 


, n — 1 ,n(zi 


( 1 . 6 ) 


- 3 - 




i + 1 


i 




Figure 1. BCFW recursion for color-ordered tree-level gluon amplitudes. 


= (pi +Pz -+-- mf = o. (1.7) 

While the trace bases, the color-flow bases and the DDM bases are conceptually simple and well- 
established in the literature, they all suffer from being non-orthogonal (unless N c — y oo), and - in an 
N c = 3 sense - non-minimal. The bases being non-orthogonal implies that the squaring of tree-level 
color structures involves (n — l)! 2 terms in the case of the trace and color-flow bases, and (n — 2)! 2 terms 
for the pure Yang-Mills specific DDM bases. Going to arbitrary order in perturbation theory, new color 
structures appear and the number of vectors (neglecting charge conjugation invariance) increases up to 
subfactorial(n) « (n!/e) in the trace bases [25]. For practical purposes this means that it is hard to treat 
the color structure of processes involving more than ~ 8 gluons plus qq -pairs exactly, and currently the 
most efficient technique for multi-parton calculations is probably to sample color structures in the color-flow 
basis by Monte Carlo techniques [24]. 

One way to cure the bad scaling involved in the squaring step is to use orthogonal bases. In the case of 
few partons, such bases, based on the decomposition of the color structure into irreducible representations, 
have been around for a while [47-52], but a general strategy for basis construction based on multiplets was 
presented only recently [25]. Below we give an overview of their key properties. 

• Multiplet bases 

These bases are based on subgrouping sets of partons and forcing the parton sets to transform under 
irreducible representations of SU(IV C ). By applying the same subgrouping procedure to all basis 
vectors, orthogonal basis vectors are obtained. 

The decomposition into these bases (to any order in perturbation theory) may be written 

A4( 5l , ff2 ,..., 5n )~^V“^, (1.8) 

a 

where a is some (collective) index of the basis vectors describing the involved representations and 
A a are the kinematic factors, clearly not equal to the A(a) in eqs. (1.1-1. 3). 
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As the multiplet bases have no direct connection to tree-level color structure, they do not typically 
span a minimal set for tree-level color structure alone. Instead they are applicable to any order, 
and are (or can trivially be made) minimal for any finite N c , leading to a significant reduction in 
dimension for a large number of partons (see table 1). 

As the multiplet basis vectors V“ do not (generally) have any cyclic symmetry under exchange of gluon 
indices, one can not expect recursion relations as simple as for the trace basis to hold. On the other hand, 
being orthogonal, multiplet bases speed up the squaring of amplitudes very significantly. We thus expect 
to gain in the squaring step, at the expense of a more intricate color decomposition. This decomposition - 
in a recursion context - is the topic of the present paper, and we will discuss two different ways of achieving 
it. 

• Clearly one strategy is to express the basis vectors in any of the bases where the recursion is known in 
terms of the relevant multiplet basis. While this strategy - which in principle is nothing but a change 
of basis - is straight forward, the gain in computation time is unclear as it involves decomposing 
(n — 1)! or (n — 2)! color structures into an exponentially growing number of basis vectors. In theory, 
the exponential times the factorial does scale better than the factorial square involved in squaring 
amplitudes in the non-orthogonal standard bases. In reality, however, this difference only becomes 
significant for a relatively large number of gluons, (cf. table 1). Directly rewriting recursion results 
obtained in other bases may, however, still be beneficial if one color factor (for example one trace or 
one DDM color factor) can be rewritten as a linear combination of a small number of multiplet basis 
vectors, and the non-vanishing projections can be identified quickly. For tree-level gluon amplitudes, 
the best option is likely to use the smaller DDM basis eq. (1.3). 

In principle, for tree-level processes with only gluons it is also possible to use the Bern-Carrasco- 
Johansson relation [34] (the BCJ relation has been proved in both string theory [53, 54] and field 
theory [45, 55, 56] ) to reduce the color-ordered amplitudes in a smaller basis of (n — 3)! color-ordered 
amplitudes. However, this relation will introduce complicated kinematic factors which are functions 
of the external momenta [34]. Although, there has been extensions of the BCJ relation to cases 
with fermions e.g., [57, 58], there is no general expression for the BCJ relation in the presence of 
(fundamental representation) quarks. Since the long term goal of this work is to pursue recursion for 
processes involving quarks as well, we will therefore not further consider BCJ relations. 

• The second approach is to perform the recursion directly in the multiplet basis. Once we know 
the multiplet basis decompositions for amplitudes for fewer gluons and relations for decomposing the 
color structure appearing in the BCFW-recursion, we can derive a recursion relation for the kinematic 
factors A a via the BCFW recursion for the color-dressed amplitudes. 

In the following, we will demonstrate the decomposition using the above strategies. In section 2, we 
calculate the kinematic factors A a by evaluating scalar products. Section 3 provides a derivation of color- 
dressed BCFW recursion, followed by a recursion relation for the color structure of MHV gluon amplitudes 
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formulated in the multiplet basis. This finally allows us to derive the BCFW recursion for the kinematic 
factors A a for any number of gluons. In section 4, we conclude and discuss natural extensions. 


2 Evaluation by scalar products 

In this section, we derive the multiplet basis expansion by comparing the color factors in multiplet bases 
to those in DDM or trace bases. The general framework for this construction is as follows: 

• The color vectors (including powers of Tr) in the DDM or trace basis (or, in general, any spanning 
set in which the recursion is known) can be expanded in the multiplet basis vectors, 

c CT = ^a QfT V“, (2.1) 

a 

where the coefficients are given by scalar products of these two kinds of color factors 2 

a ary = (V a |c cr ). (2.2) 


The scalar product is given by summing over all external color indices, i.e., for arbitrary color 
structures c 1 and c 2 , 

( cl |c 2 ) = ( cl aia 2 ...T C 2 aia2 _„, (2.3) 

Ql, a 2, ... 

with cii = 1 ,,N C if parton i is a quark or antiquark and a* = 1,..., N 2 — 1 if parton i is a gluon. 

• Substituting the expression eq. (2.1) into the DDM decomposition, eq. (1.3), or the trace basis 
decomposition, eq. (1.1), and collecting the kinematic factors corresponding to a given basis vector 
V Q in the multiplet basis, the color-dressed amplitude can be stated as 


M(g 1 ,g 2 ,...,g n ) = 9 n 

a 




V". 


(2.4) 


Comparing the above expression to the multiplet basis expansion eq. (1.8), we can read off the 
kinematic factor multiplying the basis vector V“ 


A a = ^a ao A{a). 

<J 


(2.5) 


In principle, one can use this method to derive the multiplet basis expansions for an arbitrary number of 
external legs with arbitrary helicity configuration. We have calculated the multiplet basis decompositions 
for up to six gluons, and demonstrate the calculations in the three- and four-gluon cases below. The five- 
and six-gluon cases are treated using multiplet basis recursion in section 3. 

2 We assume the orthogonal multiplet basis to also be normalized, if it is not, this is trivially accounted for. 
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2.1 The three-gluon amplitude 

For three gluons the multiplet basis can be chosen as 


V 8s 
V 8a 


9192 93 = d919392 = ^[Tr(f ffl t 93 f 92 ) + Tr{t 93 t 92 t 93 )}, 
R 


919293 = if 919392 = ^r[Tr(t 9 H 93 t 92 ) _ 
making it orthogonal but not normalized. From these two equations, we get 

Tr (t 9 l t 93 t 92 ) = —- [ d919392 jJ9l9392j Tr(t 9 l t 92 t 93 ) = — 'yd 919392 — if9l9392j 


( 2 . 6 ) 


(2.7) 


Thus the three-gluon multiplet basis expansion is given by 

3 


1 yt Ri 


M(g 1 ,g2,g 3 )=g(-1=') ^[d 939392 + if 939392 ]A(l, 3,2) + g I J ^[d 939392 - if 939392 ]A(l, 2,3) 

= g-L=d 939392 [A(l, 3, 2) + A(l, 2,3)] + g^Y f 919392 [A(l, 3, 2) - A(l, 2,3)]. (2.8) 


Since 7l(l, 2, 3) is antisymmetric under 1 2 we find 


^8s = 0 A 8a = -=A(1,3,2), 


(2.9) 


where we note that the first equation can be seen as a manifestation of charge conjugation invariance 
(cyclic reflection), and that the second color factor is precisely the amplitude for the (only) vector in the 
DDM basis. 


2.2 The four-gluon amplitude 

In the four-gluon case, we start from the DDM decomposition with 1 and 2 as the fixed legs 

M(g u g 2 ,g 3 , g 4 ) = g 2 ^if 9l92h if 19394 A(l, 2,3,4) + 3,2,4). (2.10) 

Using ColorMath [59] to evaluate the scalar products in eq. (2.2) this is decomposed into the multiplet 
basis 3 * V 


Mvj 

given by [25], 


V : 

i > v ■ 


8s 


„V 


8 a 


v 27 — 

91 93', 92 94’ 91 93i 92 94 ’ '3133; 92 94 5 3i 93 ; 32 34 ’ ' 3i 93', 92 34 1 ’ 3i 93', 92 34 


1 


V 


10 


+ V 


10 


1 


,V] 


o 

31 33; 32 34 


yd _ pa 

V 9i 93', 92 94 91 93; 92 94 


},(2.H) 


( 2 . 12 ) 


3 The gluon order convention in V gig3;g2g4 probably seems unnatural at this stage, but the advantages will become clear 

in section 3.3. 
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where P Q is the projection operator onto the irreducible representation a with dimension d Q [25, 50, 60-62], 


91 9 3 ; 92 9 a jy2 _ ^ 


§9 1 93 $92 9 4 


>8s 


N r 


91 93; 32 34 2Tr(N£ — 4) 
-1 


^9193*1^19294 


>8a 

91 93; 92 94 


p 2 7 _ 1CA91 *1X93 *2 i £91*2 £93*1 j 

9193; 92 94 4W u ~ u u > 


2N C T R 


l jig 193 *1^ jrn 92 94 


1 


5 * 192^294 + J_ t r(flt 94 f 2 t 92 ) 

■Tr 


_^c__2 8s 

2_/\r r 9i 93; 92 94 


W - 1 


2jv c 9193; 92 94 


plO — HfA9l*l£93*2 _ £9l*2£93*l j 

9193; 92 94 4^ ' 


plO — _ CA91 *1£93 *2 _ £91 *2£93 *1 'i 

9193; 92 94 4^ > 


pU _ _ ( X9l *l£93 *2 1 £91*2 £93 *15 

9193; 92 94 4^ ° u * 


5 * 1925*294 + i_ tr (fit94f2 t 92) 
T R 

5 * 1925*294 _ J_tr(f lt 94 t *2£92) 
5 * 1925*294 _ A.tr(t* n 9 H i2 t 92 ) 

t r 


_ p8a 

2 9193:9294 


_ p8a 

2 9193:9294 


+ 2 8s + 1 D 1 


25y c 9193:9294 


2jv c 9193:9294 


(2.13) 


and the general expressions for the dimensions are 1 

4 = N c - 1’ 4o = ^(iV c 4 -51V c 2 + 4), d 2 7 = ^c 2 (^c+2iV c -3) d 0 = ^N 2 (N 2 -2N c -3). (2.14) 


Expressed in this basis (which is also electronically attached as an online resource) the amplitude can be 
stated 


A 4 (c/1,32, 23 , 34 ) = g 2 A- V, 


(2.15) 


where A is the kinematic factor, 

A = N c x {-24(1, 2,3,4), -V(N C + 1) (N c - 1)A(1, 2, 3,4), -y/(N c + 1) (N c - 1) (2.16) 

[-4(1, 2, 3,4) + 2A(1,3, 2,4)] , aJ{N c + 3) (N C -1)A(1, 2, 3,4), 0, - y/(N c + 1) (N c - 3)i(l, 2,3, 4)}. 

Note that in the above discussion, we did not specify the helicity of the external legs. When we want to 
consider the kinematic factor for a particular helicity configuration, e.g., l~,2 + ,3 + ,4 _ we just substitute 
the corresponding form of the color-ordered amplitudes into eq. (2.17). 

independently of N c we refer to the adjoint representation as the octet representation, and similarly we label other 
representations by their N c = 3 dimension, although clearly the dimension depends on N c . 






















A few remarks on the basis are in place. First we note that for N c = 3 the last basis vector vanishes as 
it corresponds to a multiplet which only appears for N c > 4. For four gluons this reduction in dimension 
due to small N c is rather unimportant, but for large n, the difference becomes significant, (cf. table 1). 

Then we note that, due to charge conjugation invariance (which at tree-level manifests itself as cyclic 
reflection in trace bases), the decuplet and anti-decuplet in eq. (2.11) must multiply the same amplitude, 
which - at tree-level - vanishes. Charge conjugation invariance is also the reason why the octet vectors 
corresponding to contracting one antisymmetric structure constant with one symmetric structure constant 
vanish. This means that the four-gluon basis vectors can be expressed in terms of projectors only. Note, 
however, this a special feature of the four-gluon basis, it is not generally possible (even for n even), as there 
are different ways of building up the same multiplet, see for example [25]. 

Finally we point out that although the multiplet decomposed result eq. (2.17) may look somewhat 
complicated, it is now in an excellent form for squaring. For given helicities and external momenta, the 
color-ordered amplitudes in eq. (2.17) are just (complex) numbers. To get the full amplitude square, we 
thus just have to square the coefficients in eq. (2.17) and add them up. 

In this particular case of four gluons only, not much is gained by this rewriting, as the scalar product 
matrix for the DDM basis anyway only involves 2x2 terms. However, for larger bases, where the 2x2 
scalar product matrix is replaced by a matrix of dimension (n — 2)! x (n — 2)! (or (n — 1)! x (n — 1)! for 
trace bases), avoiding the factorial square is clearly desirable. 

Unfortunately, with the naive way of calculating scalar products utilized here, involving direct eval¬ 
uation of (n — 2)! x (the number of multiplet basis vectors) entries, what is gained in the squaring step 
for multiplet bases, may be lost in the step of scalar product decomposition. We do note, however, that a 
more clever procedure for evaluating scalar products, based on the birdtrack method and Wigner 3 j and 
6 j coefficients [26, 62] possibly could change this conclusion. As it is unclear if the scalar product method 
is beneficial, the remainder of this paper instead focuses on deriving recursion relations directly in the 
multiplet basis. 

3 Recursion in multiplet bases 

In this section, we present an on-shell recursion approach for the kinematic factor A a in the multiplet 
basis expansion eq. (1.8). We will show that once we know the recursion relation between the color factors 
in the multiplet bases for the n-gluon and (n — l)-gluon amplitudes in addition to the BCFW recursion 
for color-dressed amplitudes, we can derive a recursion relation for the kinematic factors A Q for the MHV 
helicity configuration. The main idea is: 

• We use BCFW recursion to rewrite the color-dressed amplitude in terms of products of on-shell lower 
point color-dressed amplitudes for all MHV channels. 

• For a given channel in the BCFW recursion, the color factors in the BCFW expression of the n-gluon 
amplitude can be constructed by vectors in the (n— l)-gluon multiplet basis contracted with structure 
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constants, while the corresponding kinematic factor for the BCFW expression is obtained from the 
(re — l)-gluon on-shell MHV amplitude and the three-gluon MHV amplitude. 

• We derive a recursion relation for the color structure between the n- and (re — l)-gluon multiplet 
bases. Using this recursion relation we can express the color factors in the re-gluon multiplet basis, 
and collecting the kinematic factors corresponding to the same multiplet basis vector, we obtain a 
recursion relation for the n-gluon kinematic factor A a . 

In the following, we first present a review of BCFW recursion for color-dressed amplitudes and then discuss 
the MHV configuration. 


3.1 BCFW recursion for color-dressed amplitudes 

Now let us review the BCFW recursion for color-dressed gluon amplitudes at tree-level [63]. We consider an 
re-gluon color-dressed tree amplitude M(gi, g%,..., g n ), where gt is used to denote a gluon with momentum 
(counted outwards), helicity and color. If we shift the momenta of the gluons g\ and g n with a complex 
four-vector q obeying 

q-pi = q-Pn = q 2 = 0, (3.1) 

the shifted momenta, 


pi(z) =pi - zq, p n (z)=p n + zq, 


(3.2) 


remain on-shell. With this shift, the color-dressed amplitude A4(z) becomes a rational function of the 
complex variable z. The desired amplitude is just Ad (0). To solve for A4(0), we use Cauchy’s theorem 


/ 


dz 


M(z) 


finite poles 



2:^00 


(3.3) 


The integrals around the finite poles are given by their residues, thus we have 

M(0) = -J2 Res^^^ + B, 

ZiAO 


(3.4) 


where B comes from the contour integral at infinity. In the study of BCFW recursion, the boundary 
behavior when z —> oo is important and has been investigated systematically in [64]. For gluon amplitudes, 
we can always choose a shift such that B = 0. The residues of the finite poles can be obtained by considering 
the factorization behavior, with which the amplitude is factorized into two on-shell sub-amplitudes when an 
internal line goes on-shell. The nontrivial contributions for the z-poles are those with the two shifted legs 
in two different sub-amplitudes. For the shift in eq. (3.2), we let the gluon g\ be in the left sub-amplitude 
and the gluon g n be in the right sub-amplitude. If we divide the other (re — 2) gluons into the left set T 
and the right set J , as in figure 2 , the position of the pole corresponding to this division can be found by 
solving 
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1 


n 


1 



J 


Figure 2. The amplitude is divided into a left set X U 1 and a right set J U n. Note that here, as opposed to in 
figure 1, no ordering is inferred among the partons on the left and right side. 


P X (z) = Ipi( z ) + Y Pk ) =0, 

V fcex / 


(3.5) 


giving the poles 


z X = 


Pi + E Pk 
kei 


(3.6) 


2? • ( Pi + E Pk 
\ hex 

The n-gluon color-dressed amplitude is factorized into two lower point on-shell amplitudes 


M(z) Y M (gi(z x ),{g k \kex},-Pz ,hx {zxj) g J M (P hx (z x ),{gi \iej},gn(z X j) , (3.7) 

ix,hx ”xv z l) 


where i x and h x are used to denote the color and helicity indices for the internal line. Then at z x 
M{ z Z)z ^ Y M (UzzUgk \kex},-Pi z ’ hx (z x )) ^ - M (p?’~ hx 


( z x ),{gi \iej},dn(z x ) ) . (3.8) 


ix,hj 


' z x p£(z x ) 

Considering the position of the pole in eq. (3.6), we have 

M (gi(z x ),{g k | k £ X },-P x ’ hl (z x ) \ M (P^~ hl (z x ),{gi \i eJ },g n (z x )\ x 


M(z) 


■E*- 

ix,hx 


Pi 


Z- z x 


,(3.9) 


where the denominator is the squared sum of the momenta of the external legs in the left set. Thus the 
residue for the division X, J is 

M(z) 

■ = - JVI [yn z u>Wk \k&xh~ r x j- 


Res 

Z-*Zx z 


Y M (^1 i z x),{gk Uex}. -P x ’ hx {z x )^ -^M (Pp' hx {z x ), {gi \i&j},g n (z x )y 


ix,hx 
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(3.10) 


The numerator here is a product of two on-shell sub-amplitudes with the momentum of gluon g\ and 
g n shifted. To get the full amplitude, we should sum over all possible residues at all finite poles. The 
amplitude is then given by the following BCFW recursion relation 


M ( 01 , 52 , • • ■;9n) 

= X] M (si ( z i),{9k \k&z},-Pz' hT ( z i)^-j52 M ( P£~ hl ( z z),{gi | i £ j},g n (zz)y (3.11) 

Z ixXx 

3.2 Kinematic recursion 


When we calculate an (n > 3)-gluon color-dressed amplitude for a given helicity configuration, the con¬ 
figurations with all helicities positive and all helicities except one positive have to vanish [6]. The first 
nontrivial configuration is the MHV configuration with two negative helicity gluons. 

It is convenient to use the spinor helicity formalism [1-5] to study amplitudes. In the spinor helicity 
formalism, one expresses external momenta p% by double spinors (Ai) a (Aj)d, where (A i) a and (Aj)d are 
two-dimensional Weyl spinors. The polarization vectors are explicitly expressed as £^ + r^j — and 


~ where the spinor products are defined by (A | p) = e ba X a pb and 

[Mb 1 ] 

p and p to denote the reference spinors. The matrices e ab and e ab are given by 


A | /i 


= e“ fe A, 


^ab ^ab 



aPj,, using 


(3.12) 


In the spinor helicity formalism, all amplitudes are expressed by spinor products. The color-ordered MHV 
amplitude is given by the famous formula 

v4 


H(l+,2 + ,...,r,...,j-,...,n+) =* 


(i I 3) 


(1 | 2) (2 | 3)... (n- 1 | n) (n \ 1) : 


(3.13) 


which was conjectured in [28] and proven in [6]. When we consider the MHV amplitude where all helicities 
are flipped, we just (up to a factor (—l) n ) replace ( | ) by [ | j 


Afl-,2" 




.,j+,...,n-) 


[i I il 4 


(3.14) 


[1 | 2] [2 | 3] ... [n - 1 | n] [n \ 1]' 

Let us now consider the color-dressed MHV amplitude with g\ and g n as negative helicity gluons. We 
can conveniently shift the momenta of g± and g n in the spinor helicity formalism 


Ai — > Ai, Ai — > Ai — z \ n , 

An A n , A n —)• A n + ^Ai- (3.15) 

This is nothing but the spinor expression of the (l,n)-shift defined by eq. (3.2), as p\ ~ A 1 A 1 , p n ~ X n X n , 
and q ~ AiA n . The constraint equations, eqs. (3.1), forcing the momenta in eq. (3.2) to remain on-shell, 
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Figure 3. The allowed divisions, (3, n— 1) (left) and (n — 1,3) (right), for the MHV amplitude with the gluons with 
shifted momenta 1 and n as the negative helicity legs. As is proven in appendix A, only the {n — 1,3) factorization 
channel contributes. 

are automatically satisfied in this spinor expression and the boundary contribution B vanishes, meaning 
that we can use the BCFW recursion for the color-dressed amplitude eq. (3.11) without any boundary 
correction. We also note that with this shift, g\ and g n must be in opposite sub-amplitudes. For MHV 
amplitudes, there are only two types of divisions, sketched in figure 3, which possibly could contribute 

• The (3, n — 1) divisions with 3-gluon MHV amplitudes as left sub-amplitudes and (n — l)-gluon MHV 
amplitudes as right sub-amplitudes. 

• The (n — 1, 3) divisions with (n — l)-gluon MHV amplitudes as left sub-amplitudes and 3-gluon MHV 
amplitudes as right sub-amplitudes. 

The other divisions always contain sub-amplitudes (for more than three gluons) with less than two negative 
helicity gluons and have to vanish. In fact, as is proven in appendix A, the (3, n — 1) division also vanishes. 
Thus the full color-dressed amplitude can be stated 


M (g ] ,g£,-- 

n— 1 

> 9n—1 ’ 9n ) 



S M > ■ 

1=2 

Q + p l i~ „+ 

• • ) iJi-l-i 1 i t n > yi+v ■ ■ 

1) 7 

Jr i,n 

(3.16) 


where i* denotes the (implicitly summed over) color index connecting the sub-amplitudes. The right 
sub-amplitude is given by 


M ( -P& 


i 9i ) 9n ) 9 


WTr 




MHV 


' — P + 

> i,ni 


n 


= 9 


VT^ 




i I ~Pi-, 


-Pi 


n 


,(3.17) 


n i\ 


The left sub-amplitude is given by the (n — l)-gluon multiplet basis expansion 


M [9l , ■ ■ ■, gtv P !'n > 9i+V • ■ • > 9n- 1, 
g n-3 v“ (l-, ...,(*- 1 ) + ,P- n , (i + 1) + , . . • , (n - 1) + ) , 


(3.18) 
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where, at this point, we make no statement about what gluons are counted as incoming and outgoing in 
the multiplet bases. The recursion relation between the n-gluon color factors and the (n — l)-gluon color 
factors can be written as 

u**- = E <*«)*. (*•«» 

0 

The matrices Tj describe the effect of emitting gluon n from gluon i (from the vector a in the (n — l)-gluon 
basis), decomposed into the n-gluon basis. We will refer to these matrices as the radiation matrices for 
(n — 1) — > n gluons, and in section 3.3, we will show how to calculate them efficiently. Inserting eq. (3.18) 
and eq. (3.19) into the BCFW expression for the color-dressed amplitude, eq. (3.16), and collecting the 
kinematic factor corresponding to 'Vg 1 ...g n , we obtain the recursion relation for the kinematic factor in the 
MHV configuration 

(l-,2+,3+. 

n —1 

= E E ( T <w x 

i=2 a 

where Sj, n = ('Pi + p n ) 2 = %Pi • Pn- 

The above recursion relation for the kinematic factors expresses the kinematic factor of the n-gluon 
MHV amplitude in terms of the (n — l)-gluon MHV amplitude and the three-gluon MHV amplitude. 
To calculate the n-gluon kinematic factor for the MHV configuration using eq. (3.21), we thus use the 
kinematic factors in the multiplet basis expansion of the (n — l)-gluon MHV amplitude and the matrices 
Tj (which will be derived in the next section) as input. 

3.3 Color structure recursion 

Before stating expressions in the multiplet bases we need to fix our conventions. The vector space of interest 
is the overall singlet space for the involved (incoming plus outgoing) partons. Clearly the basis vectors for 
this space can be chosen in many different ways. The prescription detailed in [25] constructs vectors by 
first constructing gluon projection operators projecting on irreducible representations for [n/2j —» [n/2j 
gluons. Following this, basis vectors for processes with up to n gluons can be constructed. (The extension 
to processes involving quarks is achieved by grouping the quarks and antiquarks to qq- pairs and noting 
that each pair transforms either as a singlet or as an octet.) For a process with n gluons, the gluons are 
divided into [n/2] “incoming” gluons and [n/2j “outgoing” gluons, such that there are either equally many 
outgoing and incoming gluons or one more incoming gluon. For the full set of gluons to transform under 
a singlet, the overall representation under which the “incoming” gluons transform must match the overall 
representation under which the “outgoing” gluons transform. The total dimension of the vector space is 
thus given by the number of ways of combining matching “incoming” and outgoing representations. 


,n 


(3.20) 


(1~, • • •, (i - 1) , P~ n , (i + 1) ,..., (n - 1) + ) -(-»)- 

Si.n 


1 

i | Pi,n 

3 

VTr 

Pi,n | Tl 

[n | i\ 
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The gluons on either side are then subgrouped such that the first two gluons transform under rep¬ 
resentation ai, the first three gluons transform under representation 02 , etc. The set of representations 
is collectively referred to as a , and the “incoming” gluons are taken to be 51 , 53 , 52 [|]-i) whereas the 

“outgoing” are given even numbers 52 , 54 ,..., 52 • Using these conventions, and letting single lines denote 
the adjoint representation and double lines denote arbitrary representations, the orthonormal basis vectors 
are 



Here the vacuum bubbles in the denominator are Wigner 3 j coefficients. They can be normalized to one, 
inducing a normalization for the generalized vertices connecting the representations. 

Letting A a { 1,2,... ,n) denote the amplitude (for convenience we keep the gluon arguments in A a in 
this order) the decomposition into these bases may thus be written 




= 9 


n—2 


Vv“ 

/ > 91 93 •■■92r§1—09294 •■■92LJJ 




. ,n). 


(3.23) 


A remark on the implication of charge conjugation is in place. As gluons transform under the charge 
conjugation invariant adjoint representation, any overall gluon amplitude must respect this symmetry. This 
is manifest in the DDM bases in the sense that each spanning color structure obeys this symmetry, but it 
is not manifest in the trace bases and the multiplet bases. For tree-level trace bases charge conjugation 
invariance instead shows up as cyclic reflection. For multiplet bases, charge conjugation invariance displays 
itself by the amplitude for a (non-invariant) basis vector and its conjugate being equal up to a sign - and 
by the vanishing of amplitudes for which all involved representations are invariant, but the invariance is 
spoiled by the generalized vertices. 

For many gluons almost all of the basis vectors contain at least one representation which is not charge 
conjugation invariant, meaning that almost every basis vector must occur with its conjugate. Using these 
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linear combinations as basis vectors thus reduces the dimension of the vector space by approximately a 
factor two. 

For the explicit calculations for four, five and six gluons, we have used conjugation invariant bases. 
However, for comparison, the dimensions of the vector spaces, are - as for the trace basis case - stated 
without this symmetry in table 1. This also has the advantage that the vector space dimension for n gluons 
is approximately equal to the dimension for n — n q gluons and n q (/(/-pairs, see [26]. 

With the above basis conventions the radiation matrices, eq. (3.19), are given by 


V" 

91 93 ••• ffp |- n —1 i 


i; 92 g4 • 


• ^2 L ri-l j 


9i—*H 


ifgikgn = ^ (Tj) 
P 


)/ 3 a g i 33 ■ 


•92rf]No9 2 94-92|ftJ ' 


giving for the amplitudes, eq. (3.21), 


71—1 

^ (l-,2+,3+,...,n-)=^E( T ‘W 


( 9+ 1)+ (V_|_li+ ( n + 1 1 


i | Pi,n 

3 

^1 , - , . . • , {l 1) L) j ^ 

Pi,n | ^ 

[n | i ] 


(3.24) 


(3.25) 


As can be seen in eq. (3.24), the color structure of the recursion relation in the multiplet bases is given by 
inserting one gluon to the (n — l)-gluon basis vectors, 


[ v“ 

l 3133 ■■■3n[-w-l- [ 1 ;g2g4-S^rt-lj 


9i—*ii 


i f9iU9n 


(3.26) 


For example, denoting the five-gluon basis vectors V“ l933s;g294 , if we radiate another gluon g$, we can 
attach it to any of the gluons g i, (/ 2 , 53 , 54 and 55 . The corresponding color factors then become 


V" i f9im 

*33 35; 32 34 J 


V* 


j f 32*36 ; Y° 


i J33*36 ; 


V“ 

31 33 35 i 32 


.if94*96 anc J Y“ 


31 S3 35; * 34 J 1 31 *35; 32 34 J ’ 31 33 35; 32* J 31 33 *; 32 34 

systematic evaluation of such color structure in the larger basis is the topic of the present section. 


ifgsm. The 


One way of evaluating the radiation matrices is to simply calculate scalar products between the left 
hand side color structure of eq. (3.24) and the vectors in the larger basis. This is equivalent to the method 
of section 2. However, as most of these scalar products vanish - for reasons that will become clear later in 


this section - such a strategy would be unnecessarily expensive. 


Instead we here present a more elegant way of evaluating the weights for the n-gluon basis vectors 
using group theory and the birdtrack notation [62] . By applying group theoretical relations to the left hand 
side of eq. (3.24) it can be cast into the form of the right hand side, with the radiation matrix elements 
expressed in terms of group theoretical weights, the Wigner coefficients. 


A method for evaluating scalar products between Feynman diagrams and multiplet basis vectors is 
explored in [26]. The same techniques are applicable for the decomposition of an (n — l)-gluon multiplet 
basis vector which has radiated an nth gluon, into n-gluon multiplet basis vectors. For this decomposition, 


three group theoretical relations are required, the completeness relation for tensor products, the color 
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structure of a vertex correction and the relation between the ordering of the representations of a vertex. 
In birdtrack notation the completeness relation reads 



In the tensor product /i <S> v above, there can be more than one instance of a particular representation. 
In this case all instances have to be summed over, for example in A ® A, where A denotes the adjoint 
representation (not to be confused with the amplitude), there are two “octets”. 

The second relation is a special case of the Wigner-Eckart theorem, the color structure of a vertex 
correction can be written as [26, 62] 


c 




(3.28) 


The above equation is sometimes stated without the sum. Indeed the sum is only needed if there is more 
than one instance of 7 in the tensor product a® (3. In this case, the a/3 7 -vertex may appear in more than 
one version, this occurs, for example, when a, f3 and 7 are octets, then a is if and d. 

The third relation concerns the ordering of representations in a vertex, the relation between the two 
orderings is given by [62], 




(3.29) 


where the equivalence defines Yutsis’ notation [65]. Typically this just gives a sign ±1, for example we 
have a minus sign for the antisymmetric triple-gluon vertex. 


3.3.1 Example: 4—^5 gluons 

The method of evaluating the radiation matrices with the above stated relations will first be applied to a 
4 —> 5 gluon example, and after that a general formula will be derived. Let us thus consider a four-gluon 
basis vector radiating a fifth gluon from gluon 3. In diagrammatic form, denoting the standard triple gluon 
vertex, if abc 1 with a black dot, where the indices are read in counter clockwise order, the color structure 
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becomes 


1 


2 



N ai 


where we have drawn the fifth gluon such that the gluon ordering in eq. (3.21) is respected. 
Applying the completeness relation eq. (3.27) to gluon 5 and the representation aq gives 


(3.30) 



(3.31) 


where the sum runs over representations in the adjoint representation times aq. On the right hand side 
above, gluon 1 and 3 and the representation f3\ are connected by a vertex correction. Using eq. (3.28) with 
7 —> /3±, e —> ctq and a, /3,5, ( —>■ A to remove it, the color structure is 




This equation is now of the desired form, eq. (3.24) where the radiation matrix components trivially can 
be read off by comparison to eq. (3.21). Letting (3 = (/3 i,/32 ) denote the representation set labeling the 
5-gluon basis vector, we immediately see that the representation /?2 is constrained to be @2 = aq. Thus 
most of the projections onto the 5-gluon basis vectors vanish. 

Note that eq. (3.32) has been derived without explicitly stating the representation aq. The result is 
thus generic and, knowing the Wigner coefficients and the dimensions of the representations, it can be 
used for immediately writing down the decomposition for any initial V“. As an example, if aq = 10 the 
allowed j3\ representations, i.e., those present in both A®ot\ and A® A, are 8, 10, 27 and 0 (for N c > 4). If 
(3\ = 8, there are two possible vertices connecting gluon 1 and 3 in the 5-gluon basis vector, and similarly 
if (3\ = 10 (and N c > 4) there are two vertices connecting /3q and a,\. 

For evaluating the right hand side of eq. (3.32) we use the dimensions of the representations, stated 
in eq. (2.14), and the Wigner coefficients calculated as in [26]. Ordering the allowed representations (3\ as 
(8 s , 8 a , 27, 10-f, 10^ d ,0), the Wigner 6 j coefficient in eq. (3.32) takes the values 

V2 


-1 


. 0 , 


1 


-, 0 , 


1 


V / iVF r 4(A r c 2 -l)’ ’ V N c + SN c + 2 (1V C 2 - N c ) ’ y/.W ^4 (AT2 - 1) ’ ’ yW c 2 - 3!V C + 2 (1V C 2 + N c )) 

(3.33) 
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respectively. In [26] the Wigner 3 j coefficients are normalized to one. Eq. (3.32) is valid for any nor¬ 
malization as long as it is consistently used. However, requiring all 3 j coefficient to be one implies a 
normalization if abc if cba = 1 for the triple-gluon vertex. To get the correct i/ a 6 c -normalization, eq. (3.33) 
must be therefore multiplied by a factor of 



y / 2iV c (A^| — 1 )Tr (using standard normalization of vertices). 


(3.34) 


Normalizing the Wigner 3 j coefficients to one and using the definition of the normalization constants, 
eq. (3.22), gives 

N ai 1 

-rr -5 - = (3.35) 

Combining the dimensions, the Wigner 6 j coefficients, the normalization factor eq. (3.34), and the overall 
sign in eq. (3.32) gives 



( I 2N C jN C (N C + 3) 

\ V N% — 4’ U ’ y 2(N C + 2) 


— \fN~c, 0 , 


N C (N C - 3) \ 

2(Ay - 2 ) J ' 


(3.36) 


These are the factors required to express the color structure of eq. (3.30) in terms of the five-gluon basis 
vectors. For bases which are not charge conjugation invariant, they would be the entries in the radiation 
matrix T 3 , corresponding to mapping the initial vector V 10 emitting a gluon g$ from gluon g 3 onto the five- 
gluon basis vectors V 8s,1 °, V 8a,1 °, V 27,10 , V 10,10 ^, V 10,10 ^ and V 10: °. In the above case, the initial color 
structure, eq. (3.30) with oq = 10, is from a basis vector which is not charge conjugation invariant, and 
neither is any of the vectors which the color structure is projected onto. To get to the charge conjugation 
invariant vectors in this case simply requires the substitution of 10 —> 20 in the vector representation labels. 
With this change, the result, eq. (3.36), can be compared to column five of T 3 in eq. (B.3) in appendix B, 
where the radiation matrices for 4 —> 5 gluons are given expressed in the basis from eq. (B.l) (which is 
also electronically attached as an online resource). 

In general, when knowing the radiation matrices in non-charge conjugation invariant bases and con¬ 
verting to charge conjugation invariant bases, a sign might be required since the n-gluon non-invariant 
vector may come with a minus sign in the linear combination building up the charge conjugation invariant 
vector. For the same reason another minus sign may come from the (n — l)-gluon basis vector. Apart from 
a potential sign, factors compensating for vector normalizations and occurrence of both a vector and its 
conjugate on the right hand side of eq. (3.24) may be required. 

For five external gluons (and N c > 4) there are 22 charge conjugation invariant basis vectors, stated in 
eq. (B.l), but the example color structure is given by only four of them. It is worth noting that although, 
in this case, most of the basis vectors with /?2 = 10 contribute (there are only two zeros in eq. (3.36)), for 
more external gluons the constraints corresponding to /3q £ aq < 8 > A become more restrictive. This point 
will be elaborated on after the derivation of the general formula for radiation matrices. 
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3.3.2 The general case: n — 1 —> n gluons 

In general, if the nth gluon is radiated from one of the incoming gluons, the color structure of eq. (3.26) 
will be of the form 


l 

3 
5 

jy Q l a2...Q n -l-3 

i = 2j + l 

2 (j+l)+l 

2(Pr i l- 1 )+ 1 

where, if j = 0 or 1 (i.e. i = 1 or 3), otj-i is an octet. In the following steps it will be assumed that gluon 
1 is not the emitter, this special case will be addressed after the derivation. To get to the right hand side 
of eq. (3.24), the same steps as in the above example are used: First the completeness relation, eq. (3.27), 
is applied repeatedly and then vertex corrections are removed using eq. (3.28). 

We thus want to apply the completeness relations to gluon n and the representations a.j. ctj + 1 ,..., Q!|-n = i-|_ 1 , 
i.e., we insert it in the encircled positions in 



(3.37) 


jy-Ql 02 ... 





resulting in 


(3.38) 
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Here gluon 2 j + 1 has a vertex correction of a different form compared to the other gluons. We first 
treat this separately and then address all other vertex corrections. Using Yutsis’ notation, eq. (3.29), and 
eq. (3.28), the leftmost vertex correction can be written 



The remaining vertex corrections can be contracted similarly, for example 



(3.40) 


(3.41) 
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Applying the above steps and adjusting the vertex order of the last new vertex results in 


E 



“fc-i 





■Wi-i 

/a r^i-i 


Wl 


16 rvi 


N^i 


Q!2...0n-4 


-N 


Qt-1 ...Qtj-i(3j . ../3 j- n_l -| _ 1 ®i-n- 1 1 _i •••Ot-n—4 


jy® 1 ' 1 ’''^\7 — 1 • • • /^ |- 7T- — 11 _ ^ ^ |- n — 1 — 4 

i Pj Pj+ 1 2 


<*i a2 


Ql rn-ln 

I-2 - ! 



(3.42) 


where the vertex labels in the third line have been suppressed. The third line of eq. (3.42), combined with 
the last normalization factor of the second line is now of the form of a basis vector. Hence the factor in 
front of it is the radiation matrix T 2 j+\ expressed in terms of Wigner coefficients. 

We remark that the form of the radiation matrix from the example, eq. (3.32), differs from eq. (3.42). 
This is only due to the fifth gluon being drawn in an, at that point, more natural way. The two expressions 
are identical, if the expression from eq. (3.42) is written out for gluon 3 with n = 5, it can be simplified to 
become exactly the expression of eq. (3.32). 

The derived result, eq. (3.42), is for gluons emitted from the “incoming” gluons. For gluons emitted 
from the “outgoing” gluons an analogous derivation can be done, resulting in an equation similar to 
eq. (3.42). A special case occurs if the emitter is the first gluon on its side, gluon 1 for the “incoming” side 
and gluon 2 for the “outgoing” side. Compared to the radiation matrices T 3 and T 4 the matrices Ti and 
T 2 are identical up to sign differences for some entries. The difference originates from a difference in the 
leftmost (rightmost for gluon 2) vertex correction, eq. (3.40), that results in a change in vertex ordering 
which gives a possible sign difference. (Taking gluon 1 and gluon n to be the gluons with shifted momenta, 
emission from gluon 1 need not be considered, see figure 3.) 
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Concerning the group theoretical constraints on the representations we note that each representation 
to the left of / 3j and to the right of /3|-n-i-| 1 are constrained to be equivalent to representations in the set 
a, i.e., 


Pk = a k , for k = 1,... ,j — 1, 

Tl — 1 

Pk = Oik- 1, for k = r ^ 1, - • • ,n - 3. (3.43) 

The shift in the index of the second equality in eq. (3.43) is from the fact that there are n—3 representations 
for n gluons and n — 4 for n — 1 gluons, with the new representation being inserted just before aj-n-i-j 1 . 

There are also constraints on the set of representations ft coining from the completeness relations 
applied in eq. (3.39), 

72 — 1 

Pk^Oik® A, for k = j,j + 1, • • •, r— 1 - !■ (3.44) 

This constraint will become more restrictive when the emitting gluon is far from the middle of the basis 
vector, since in this case, the constraint is imposed on many representations in the basis vector. 

Using the constraints eq. (3.43) and eq. (3.44), the number of possibly non-zero elements in the radi¬ 
ation matrix columns can be counted (more are zero due to generalized vertices giving Wigner coefficients 
which are zero). The result, averaged over all possible emitters 2,.., n — 1 and all possible initial basis vec¬ 
tors, is shown in table 1 along with the maximal number of possible P for any a. Both the maximal number 
of possible P for any a and the average over all a and all emitters are overestimates, as symmetries of the 
Wigner coefficients will force some of them to vanish, depending on the choice of vertices. In addition, 
there are Wigner coefficients vanishing due to the invariance condition of tensors under the group, see [62] 
for the invariance condition in birdtrack notation. For the calculated radiation matrices, the reductions 
due to vanishing Wigner coefficients changes the averages for N c > n from 6.9 to 3.9 and 14.4 to 8.9 for 
the n = 5 and n = 6 cases, respectively. For the gluons only case there is the further reduction due to 
charge conjugation invariance. For comparison, table 1 also shows the dimensions of the (all order) vector 
space, both for QCD and in the limit N c —> oo, and the number of vectors in the spanning sets for the 
“trace bases” and “DDM bases”. 

We note that although the average (and maximal) number of terms does increase with the number of 
gluons, compared to the increase in the dimension of the vector space, this growth is very mild, meaning 
that a smaller and smaller fraction of all basis vectors contribute. Instead of having to treat the square 
of the number of basis vectors in the squaring of the amplitude, we thus only need to treat the number 
of basis vectors, times the number of contributing emitters times the average number of terms in deriving 
the amplitude; recall that the squaring of the color structure itself is quick in the orthogonal bases. For 
example, using the trace basis for 10 gluons requires (362 880) 2 ~ 10 11 terms in the squaring step, with 
the (gluon specific) DDM basis this can be reduced to ~ 10 9 , whereas using the multiplet basis would 
require up to 19 208 x 8 x 106 ~ 10' terms for the identification of new basis vectors. For more gluons the 
difference is even larger. As our long term goal is to include processes with an arbitrary number of quarks, 
we view the comparison between the non charge conjugation invariant multiplet basis and the trace basis 
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n 

5 

6 

7 

8 

9 

10 

11 

12 

Avg 

QCD 

6.0 

10.8 

17.5 

32.5 

54.6 

106 

185 

268 


N c > n 

6.9 

14.4 

24.6 

57.9 

109 

299 

593 

1 775 

Max 

QCD 

8 

33 

33 

178 

178 

962 

962 

5 220 


N c > n 

9 

44 

44 

400 

400 

4 006 

4 006 

41 256 

Vectors 

QCD 

32 

145 

702 

3 598 

19 280 

107 160 

614 000 

3 609 760 

(all orders) 

N c > n 

44 

265 

1 854 

14 833 

133 496 

1 334 961 

14 684 570 

176 214 841 

Trace 

any N c 

24 

120 

720 

5 040 

40 320 

362 880 

3 628 800 

39 916 800 

DDM 

any N c 

6 

24 

120 

720 

5 040 

40 320 

362 880 

3 628 800 


Table 1. The average (taken over all initial vectors and the emitters 2, ...,n — 1) and maximal number of non-zero 
elements in the columns of the radiation matrices for (n — 1) —» n gluons. The stated numbers are overestimates since 
they assume that no Wigner 6 j coefficient involving admissible representations vanishes. For comparison we also 
show the total dimension of the all order vector space for N c > n and for N c = 3 (without accounting for reduction 
due to charge conjugation invariance) and the number of spanning vectors in the tree-level trace bases and the DDM 
bases. 


as the most relevant comparison. From table 1 we thus conclude that the overall treatment of the color 
structure can be sped up significantly by the usage of multiplet bases. 


3.4 The five-gluon amplitude 


Utilizing that only the (n — 1,3) division contributes, the BCFW recursion expression, eq. (3.16), for 
the five-gluon color-dressed MHV amplitude is given by the sum of the diagrams in figure 4. As seen in 
eq. (2.17), the four-gluon color-dressed (MHV) sub-amplitudes can be expressed in the multiplet basis. 

The contractions of basis vectors in the four-gluon sub-amplitudes and the structure constants in the 
three-gluon sub-amplitudes can be expanded in the five-gluon multiplet basis using the radiation matrices 
Tj for emitting the gluon <75, illustrated in figure 4. Collecting the kinematic factors corresponding to a 
given five-gluon basis vector, we can write down the recursion relation eq. (3.25) for the kinematic factor 
for five gluons 


A^(l-,2 + ,3 + ,4+,5-) 


^(T 2 ) /3q A“(T-,H 2 - )3 +,4 + ) 

a 

+ E( T S W' 4 “(V' 2+ .bT50 + ) 

a 

+ E( T ^a^(l“ ) 2 + ,3+,P 4 7 5 ) 
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1 1 

2 | 

-P2 ,5 


■S25 VTr 

2 | 

5 

5 | -P'2,5 

1 1 


3 I 

-P 3 , 5 

3 

S35 VTr 
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5 — T*3,5 

1 1 


4 | 
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3 

S45 VTr 
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5 | ~Pa ,5 


(3.45) 
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Figure 4. The relevant splittings for the BCFW recursion in the five-gluon MHV case (above) and the corresponding 
color structure for the radiation matrices (below). Note that the orders of the vertices are drawn to be 

consistent with eq. (3.25). 

Here the matrices T, in the first, second and third term of eq. (3.45) are the radiation matrices corre¬ 
sponding to the first, second and the third diagram in figure 4. They are calculated as in section 3.3 
and explicit matrices are stated in appendix B. The four-gluon kinematic factors A° 3 + , 4 + ^, 

A a (T“, 2 +, P 3 > 4+) and A a (l", 2 +, 3+ P 4 are given by replacing the gluons 2, 3 and 4 in eq. (2.17) by 
P 25 , P 35 and P± 5 respectively, and inserting the explicit expressions for the color-ordered MHV amplitudes 
eq. (3.13). Using the Mathematica package S@M [ 66 ], one can compute the kinematic factors. Clearly, in 
the five-gluon case, the only non-vanishing helicity configurations are the MHV and the MHV making our 
result applicable to all relevant cases. 

As shown in section 2 , one can always express the kinematic factor in the multiplet basis expansion 
in terms of color-ordered amplitudes in the KK basis. To achieve this in the recursion approach, we recall 
that the (n — l)-gluon kinematic factor A a { 1 , ...,n) can be written in terms of (n — l)-gluon color-ordered 
amplitudes A(a) (in the five-point case, the (n — l)-gluon result is given by eq. (2.17)). Since the (n — 1)- 
gluon color-ordered amplitudes satisfy the KK relation eq. (1.4), we can apply KK relations to rewrite the 
(n — l)-gluon kinematic factor in the multiplet basis in terms of (to — l)-gluon color-ordered amplitudes 
of the form A( 1 _ ,..., TVJ. After substituting KK expressions for the (to — l)-gluon kinematic factors - 
as well as the particular form of the MHV amplitudes eq. (3.13) - into the n-gluon recursion expression 
(eq. (3.45) for five gluons), we obtain the n-gluon kinematic factor expressed in terms of combinations of 
n-gluon color-ordered MHV amplitudes of the form A( 1 _ ,..., n _ ). Putting all this together, we arrive at 
the five-gluon result, explicitly stated in appendix B, and electronically attached as an online resource. 
Clearly, in the five-gluon case, the only non-vanishing helicity configurations are the MHV and the MHV 
configuration, making our result applicable to all relevant cases. 
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3.5 The six-gluon amplitude 

The six-gluon MHV amplitudes have been calculated analogously using the electronically attached charge 
conjugation invariant multiplet basis (see online resource). Concerning the basis, we remark that the 
dimension of the all order vector space is reduced from 265 to 140 when keeping only conjugation invariant 
linear combinations of vectors. Specializing to N c = 3 brings down the dimension further to 75. 

We also note that although - expressed in the KK basis - there could naively be up to 4! spinor terms 
multiplying each basis vector, on average only 8.5 contribute. By the scalar product approach in section 2, 
we find that the final expression for the MHV amplitude is also valid when we replace the KK basis for the 
MHV configuration by the basis for arbitrary helicity configurations. Thus we need not specify the helicity 
information in the final six gluon result. The six-gluon basis, the radiation matrices and the resulting 
amplitudes are electronically attached as online resources. 

4 Conclusion and outlook 

We have shown how BCFW recursion can be used for relating higher point tree-level MHV gluon amplitudes 
to results for fewer external legs. To achieve this recursive decomposition we have utilized two different 
strategies. 

One option is to straightforwardly evaluate scalar products of color factors in the multiplet bases with 
those in the DDM decomposition (or in principle any other basis where recursion relations are known). 
While this strategy benefits from being conceptually simple, it will be competitive for multi-particle pro¬ 
cesses only if a rapidly decreasing fraction of such scalar products are non-zero, and if the contributing 
scalar products can be identified and evaluated quickly. 

Therefore we have shown how to derive n-gluon MHV amplitudes directly in the multiplet bases. 
This requires the calculation of “radiation matrices”, describing the effect of radiating one gluon from an 
(n — l)-gluon basis vector, decomposed into the n-basis vectors. We have shown how to efficiently calculate 
these matrices using birdtrack techniques, and argued that the overall treatment of color structure can be 
sped up significantly using multiplet bases. 

While we do believe that the present paper is an important step in the direction of achieving efficient 
multi-particle amplitude calculations in multiplet bases, quite some work remains before this becomes 
reality. First of all, the results should be extended beyond MHV, to processes with quarks, and preferably 
beyond leading order. Secondly, it remains to efficiently implement multiplet bases and calculation of the 
radiation matrices. 

Finally, we remark that we have studied recursion using one particular form of multiplet bases, corre¬ 
sponding to one particular subgrouping of partons. It appears quite likely that even more efficient choices 
of multiplet bases can be made, such that the radiation matrices are even more sparse and the recursion 
can be achieved even quicker. 
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A Vanishing of the (3,n-l) division 


Here we show the vanishing of the (3,n — 1) division in the BCFW recursion expression of the MHV 
amplitude M(gi ,g^, ■ ■ ■, < 7 n) for the (l,n)-shift, eq. (3.15). The proof here is similar to the proof used 
for color-ordered amplitudes which can be found, for example, in [43]. For the (3,n — 1) division, the 
left sub-amplitudes are three-gluon MHV amplitudes, while the right sub-amplitudes are (n — l)-gluon 
MHV amplitudes. Using the fact that the sub-amplitudes for fewer gluons always can be expressed as 
linear combinations of color-ordered amplitudes in the DDM basis, we find that one contribution to the 
kinematic factor for this division has the form 


(A.l) 
. (Since there 



i\ - Pi,i 

3 /- ,_\ 3 
i . ( P M 1 n ) 

1 | i 

1 — Pl,i 



~PlA 


= ±i 


Pi 


1,2 


where l and k are two arbitrary unshifted gluons in the right set, and 
are three —Pi j ’s in the numerator and one in the denominator, the sign (± 1 ) does not appear in the 
final result.) In the above expression we have divided out a factor i\ — Pi,, in the MHV amplitude and a 
factor (^Ppj | in the MHV amplitude. The position of the pole for the above expression is 

[IN] 


\n\i\ 


si,i(z) = ([1K] - ^ [n|i]) (*|1) = 0 => zij = 

The numerator of eq. (A.l) then reads (up to a factor i) 

( i | Pi,i (Pi,i I = ([* I 1] (1 I n) ~ zi,i N I n] (1 I n}) 3 

= ([i | 1] - zi,i [i | n ]) 3 (1 | n ) 3 
= 0 3 (1 | n ) 3 . 

Also in the denominator, several factors vanish, 
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(A.5) 


Thus, after dividing out a factor ([1 | i\ — z\^ [n | i\) 2 in both numerator and denominator, the expression 
eq. (A.l) is proportional to 


[1 | i] — zi t i [n | i] = 0. (A.6) 

Therefore, the (3,n — 1) division always gives a vanishing contribution. 

B Five-gluon multiplet basis, radiation matrices and MHV amplitudes 

We use the charge conjugation invariant orthonormal five-gluon multiplet basis given below. As remarked in 
section 3.3 we need, apart from the representation labels ctj, also a label distinguishing the various vertices 
from each other. For this reason the vectors carry labels oq, «2> a4 where, ai contains additional information 
about the vertex if needed. Also, since our basis is charge conjugation invariant, the representations 10 
and 10 only appear together, referred to as 20. The five-gluon basis is 
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The radiation matrices for adding the gluon g§ to the four-gluon multiplet basis from eq. (2.11) 
expressed in the five-gluon basis in eq. (B.l) are given by 
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Using the above radiation matrices (which can also be found among the electronic attachments available 
as online resource) in eq. (3.45) we arrive at the five-gluon multiplet basis result 
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(B.5) 


where a = 1 ,..., 22, corresponding to the 22 vectors in the five-gluon basis. Here /C, C and A f are defined 
by 
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